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We show that the V2 measure of elliptic flow saturates early in the evolution 
before the hadronization transition to a rather large value ~ 0.05 — 0.15 as 
ag varies from 2-10 mb and thus is a sensitive probe of the dynamics in the 
plasma phase. 
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In nucleus-nucleus collisions at the Relativistic Heavy Ion collider (RHIC), a system of 
deconfined quarks and gluons is expected to be produced for about 3 — 10 fm/c |1[]. Various 
signatures have been proposed to verify its existence 0. Since the volume and lifetime 
of this matter may be much larger and longer than those given by the confinement scale 
~ 1/ Kqcdi it is expected that collective motion of these deconfined partons should arise 
and have observable consequences that provide information on the dynamics in quark gluon 
plasma. 

The collective motion of particles in heavy ion collisions can be studied via directed 
and elliptic flows. The directed flow, which measures the collective motion of particles in 
the reaction plane, has been studied extensively at BEVALAC, SIS, and AGS energies and 
found to be appreciable 0-^. It becomes small at higher energies due to the large beam 
rapidity. In this case, it is more suitable to study the elliptic flow [^|[, which measures 
the azimuthal asymmetry of the transverse flow pattern. For heavy ion collisions with fixed 
targets, such as at CERN SPS, elliptic flow is generated by hadronic final state interactions 
1^. On the other hand, at RHIC energies and beyond, copious mini-jet production leads 
to a high density of partons on a very short time scale and elliptic flow can be generated 
from interactions at the partonic level prior to hadronization. The aim of this letter is to 
calculate this effect using the parton cascade code ZPC P,p!0 . 



Initial partonic conditions can be estimated for Au-Au collisions at RHIC energies via 
the HIJING generator |]lT| for example. Here we will consider simple idealized geometries for 



a first analysis. Initially, partons are uniformly distributed in the overlapping region of two 
disks each having a transverse radius of 5 fm with their centers separated by a distance equal 
to the impact parameter. For collisions of large nuclei at not very large impact parameters, 
the multiplicity is roughly proportional to the overlapping area. HIJING estimates for 6 = 
a gluon rapidity density dNg/dy ^ 300. Again as a first idealization we assume a uniform 
rapidity distribution from —5 to +5 and take a momentum distribution given by a local 
thermal one at a temperature of 500 MeV. Furthermore, produced particles are taken to be 
formed on a hyperbola in the t-z plane, with z the beam direction and x-z plane the reaction 



plane, using a formation proper time tq = 0.2 fm. 

At RHIC energies, the initial minijet system is dominated by gluons. Their subsequent 
time evolution can be described by the parton cascade model ZPC. In this model, only 
gluon-gluon elastic scattering has been included, with a cross section 

dag 97r«^ /i , /"^\ 1 



1 



dt 2^5/ {t-nY' 

In the above, the strong interaction coupling constant a is taken to be 0.47, and /i is an 
effective screening mass responsible for regulating the divergent leading-order cross section. 
We vary /i to study the dependence of the elliptic flow on (Xg. 

The elliptic flow of partons is characterized by the second Fourier coefficient f 2 in their 
momentum distribution, i.e., 

1 dN 



vo + 2t;i cos(0) + 2t;2 cos(20) H . (1) 
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Alternatively, it can be characterized by V2 in the expansion of their transverse energy 

distribution, 

1 dE^ 
TEt)^ = ^0 + '^v'l cos(0) + 2v'2 cos(20) + • • • . (2) 

Numerically, f 2 is measured by averaging 

cos(20) = 4^ 

pI + pI 

over partons, i.e., it is the single particle average of squared transverse momentum asym- 
metry. Similarly, V2 is measured by the E^ weighted single particle average of cos(20). In 
general, the two azimuthal distributions do not have to be the same, and the two coefficients 
V2 and V2 are thus not necessarily equal |jl2|]. We will later give an example in which they 



develop different values. 

The elliptic flow arises due to multiple collisions as long as the initial transverse spatial 
distribution is azimuthally asymmetric. For flxed impact parameters, its value reflects the 
strength of interactions {(Jg)- For flxed cr^ its value depends on the spatial asymmetry as a 
function of the impact parameter (or multiplicity). 



In Fig. 1, we first show results for fixed impact parameter b = 7.5 fm. The time 
evolution of f 2 for different values of the gluon-gluon scattering cross section is shown. Note 
that all asymptotic values are reached very early, t < 2 fm/c, well before the hadronization 



transition. For interaction length Jcr/n larger than the mean free path, the numerical 
method of particle partition ||13[ has been used in order to obtain stable solutions of the 
Boltzmann equation. For the typical cross section between gluons, 3 mb, the final V2 after 
the partonic stage is around 0.1 which is quite "large" and easily detectable. Since the 
screening mass, or equivalently, the effective parton scattering cross section is uncertain we 
study its effect on the elliptical flow by using different values for the cross section, i.e., 1 
mb and 10 mb. We see that V2 is rather sensitive to the cross section. This demonstrates 
that very large dissipative corrections are at work since in the ideal fluid limit the flow 
depends only on the speed of sound that is in our case l/y/S. The important role played by 
dissipative phenomena in collective properties of plasma evolution was also demonstrated 
in [jl0| in connection with the transverse energy evolution. For these initial conditions even 
Navier-Stokes is inadequate and the full microscopic parton cascade is necessary. The strong 
dependence of the elliptic flow on the parton cross section could provide an important probe 
of the plasma dynamics if hadronization and flnal hadronic rescatterings do not modify the 
results too much. We cannot yet address this problem within the scope of ZPC and leave 
that important topic to a future study. 

In Fig. 2, we compare the time evolution of V2 and V2 for a gluon-gluon cross section 
of 10 mb. The statistical errors are denoted by bars. We see that V2 is always distinctively 
larger than V2 during the evolution. It indicates that particles moving in the x direction 
tend to have relatively larger Et than particles moving in the y direction. Even though 
flnal-state particle production may change the difference, this gives an example of different 
values for V2 and V2. The difference between V2 and V2 becomes less signiflcant when the 
cross section becomes smaller. 

To study how the elliptic flow reflects the collision geometry, we show in Fig. 3 the depen- 
dence of the elliptic flow on the impact parameter for medium range impact parameters. An 



almost linear dependence is observed, and this is similar to the impact parameter dependence 
obtained from the hydrodynamic model [Q . The linear increase of elliptic flow as a function 
of impact parameter reflects the initial spatial asymmetry, which can be characterized, e.g., 
by the width L^ and height Ly of the overlapping region, via as = {Ly — Lx)/{Ly + Lx). One 
can easily show that for medium range impact parameters, the spatial asymmetry is ap- 
proximately linearly proportional to the impact parameter, similar to the impact parameter 
dependence of the elliptical flow. However, for b/{2R) very close to 1, the spatial asymmetry 
approaches 1, but the V2 coefficient actually drops to 0. This is simply due to the fact that 
there are too few particles which can undergo scattering to develop a momentum asymmetry 
from the spatial asymmetry. In Fig. 4, the time evolution of f 2 is shown for different impact 
parameters. Although collisions at larger impact parameter give rise to a larger V2, it reaches 
the final value also sooner, indicating that it is the size of the overlapping region in the x 
direction which determines how fast the elliptical fiow is generated in the collision. 

To summarize, we have studied in this letter the elliptical fiow generated during the 
partonic evolution in heavy ion collisions at RHIC energies. For medium range impact 
parameters, an appreciable elliptical fiow has been found as a result of the initial spatial 
asymmetry and subsequent partonic rescattering. The magnitude of the elliptical fiow shows 
a strong dependence on the parton-parton scattering cross section, demonstrating the impor- 
tance of dissipative corrections to idealized hydrodynamic expectations. The finite parton 
elliptical fiow is expected to lead to an initial hadron matter with an asymmetric azimuthal 
momentum distribution. This is different from heavy ion collisions at lower energy, such 
as the SPS, where the hadronic evolution starts with an approximately isotropic azimuthal 
momentum distribution in the nuclear overlapping region. Since the spatial asymmetry is 
relatively smaller during the hadronic evolution in heavy ion collisions at RHIC, one does 
not expect that hadronic rescattering will generate appreciable elliptical flow 0] and thus 
should not signiflcantly modify the initial partonic elliptical flow. 

In this flrst study we have not yet investigated the effects due to different initial conditions 
for the partonic matter. For example, the turbulent glue initial conditions given by the 



HIJING model [|T4| may lead to a different value for the elliptical flow as obtained here 
using a thermal one. A quantitative study of effects due to both the hadronic rescattering 
and the initial parton distribution is underway by using a transport model that takes results 



from the HIJING model |]TT| as input to the ZPC parton cascade model and then includes 
the hadron cascade via the ART model ^. Such a comprehensive study will allow us to find 
quantitatively the sensitivity of elliptical flow on the parton-parton scattering cross section, 
thus making it possible to extract valuable information on the properties of the hot dense 
parton matter formed in the initial stage of ultrarelativistic heavy ion collisions. 
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FIG. 1. Time evolution of V2 coefficient for different effective parton scattering cross sections 
in Au-Au coffisions at -^s = 200 AGeV with impact parameter 7.5 fm. Filled circles are cascade 
data, and dotted lines are hyperbolic tangent fits to the data. 
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FIG. 2. Time evolution of V2 and V2 coefficients from the same reaction as in Fig. 1 for u = 10 
mb and b = 7.5 fm. Filled circles are cascade data with statistical errors, and dotted lines are 
hyperbolic tangent fits to the data. 
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FIG. 3. Impact parameter dependence of V2 and V2- Filled circles are for v'2 and open circles 
are for V2- The lines are splines to guide the eyes. 
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FIG. 4. Time evolution of V2 for different impact parameters. The cross section is taken to be 
3 mb. Filled circles are cascade data, and dotted lines are hyperbolic tangent fits to the data. 
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